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DECAY PROCESSES IN MHD TURBUT.ENCE: 
DYNAMIC ALIGNMENT AND SELECTIVE DECAY 


W. H. Matthaeus and D. Montgomery 


ABSTRACT 

Iftider some circumstances, incompressible magnet obydro dynamic 
turbtilence will evolve toward a state in which the velocity 
fields and magnetic fields are aligned or anti-aligned. We 
propose a mechanism for this effect and illustrate with nu- 
merical computations. Under some other circumsteuices , the 
energy appears to decay selectively toward a minimum energy 
state in which the kinetic energy has disappeared. It has 
not been possible so far to identify a boundary in the phase 
spcce which divides the two regimes. 

INTRODUCTION 

Turbulent, dissipative media obey a statisticeJ. mechanics 
whose central features are very far from clear. If the dis- 
sipative terms in the equations of motion are dropped, a trun- 
cated discretized representation of the remaining ideal equa- 
tions seems to exhibit behavior which is well described by 
classical statistical mechainics: time averages are well 

represented by phase space averages calculated by Gibbsian 
techniques. But as soon as the dissipation (viscous or ohmic, 
for exanple) is turned on with even a very small coefficient, 
the picture changes drastically. 

lypical dissipation amounts to a perfect or nearly perfect 
absorption in wavenuniber space, if the discretization is ac- 
complished by Fo\uier decomposition. The hi^-k Fourier coef- 
ficients are systematically drained, and phase points no 
longer spend times in phase volumes which are proportions^, to 
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those volumes. Even though this fractional depletion of the 
time spent in volumes of phase space corresponding to the ex- 
citations at high k is not sharply characterized, it may not 
be the only thing happening in phase space vdiich is non- 
standard, according to the perspectives of classical conser- 
vative statistical mechanics. There may also be preferences 
for the phase points of systems to spend excessive fractions 
of their time on exotic sets in the phase space: strange 

attractors, fracteuLs, and the like. Most discussion of these 
possibilities has occurred for model systems with very low 
numbers of degrees of freedom. Efforts to draw conclusions, 
from these low order "dynamical systems," about realistic 
turbulent systems with many degrees of freedom have not al- 
ways been very convincing. We have no candidates, for 
example, for a description of a set in the (Fourier) phase 
space that would lead to the Kolmogoroff spectrum. 

An alternative to the "mappings" game, which characterizes 
increasingly remote-looking model systems, is to continue to 
deal with realistic turbulent systems, emd to try for the 
present to give a physical characterization of the tendencies 
in their turbulent dynamics which appear to be genereil. This 
is also an inadequate characterization of their phase space 
behavior, but at least it may be regarded as physically use- 
ful information. 

An exeunple of such a physical effect which we wish to re- 
mark upon here is the dynamic alignment of velocity fields 
and magnetic fields in magnetohydrodynamic (MKD) tiirbulence. 
It was observed some time ago by Belcher and Davis ( 1971 ) » in 
conq>aring solar-wind magnetometer measurements with velocity 
fields inferred from particle energy analyzers, that it often 
occurred that solar wind turbulence near 1 A.U. found itself 
in a highly aligned (y = g) or anti-aligned (y = -g) state. 
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[Here y is the velocity field in a zero-ioamentum frame, ex- 
pressed in units of an rms Alfven speed, and g is the turbu- 
lent magnetic field. ] These alignment properties have been 
frequently observed since, though we attempt no survey of the 
literature here, and are often discussed rather inconclusively 
in terns of the properties of small-aiiq)litude Alfven waves. 
They have been seen also in a turbulence closure cedculation 
by Grappin et al (1982), and have been computed dynamically 
by Pouquet et al (1983) and by Ms.tthaeus et al (1982). An 
early important theoretical initiative on their role in the 
solar wind was due to Dobrowolny et al (1980). Our purpose 
here is to describe in print for the first time the computa- 
tions of Matthaeus et al (1982), and to suggest a qualitative 
e:q>lanation for the obsesrved alignment which differs from the 
one suggested by Pouquet et al (1983). 

THE DYSAMICAL PROBLEM 

We begin with the equations of incompressible MHD in a famil- 
iar dimensionless form: 

^ ♦ yVy = -Vp + + vT^y , (l) 

^ + yVg = g*Vy + yvS , (2) 

where p is the total pressure (mechanical plus magnetic) and 
where V»y = 0 = V*g. The constants y”^, v“^ are the magnetic 
and mechanical Reynolds numbers, respectively. We work in 
rectangular, periodic, two-dimensional geometry, with all 
variables Independent of the Z coordinate: 9/9z = 0. We 

solve Eqs. (1) and (2) numerically, using a spectral-method 
code of the Orszag type (Orszag 1971* Patterson and Orszag 
1971* Gottlieb and Orszag 1977) used previously (pyfe et «il 
1976, 1977 a,b; Matthaeus and Montgomery 1980, 1981; Shebalin 

et al 1983). 
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We solve the initieil value problem starting from random 
initial conditions. These are specified by choosing a set of 
random initial values for the Fourier coefficients in the ex- 
pressions 

y = y(fe.t) exp(ifc*jc) 

® exp(i^*3c) 

and then following the evolution of the Fourier-an5)litudes by 

solving tihe Fourier transformed version of Eqs. (1) and (2). 

Ihe degree of alignment may be measured by the ratio of 

two ideal invariants, the cross helicity, H , and the energy, 

c 

£: 

Hg ^/ygd^x 2j^y*(k,t)*B(k,t) 

^ |y(^,t) |^+|g(fe,t) 1^) 

5 

The ratio H^/E must lie within the range -1/2 ^ H^/E ^ +1/2, 

the limits indicating perfect anti-alignment and perfect 

alignment, respectively, and E are both conserved by 

Eqs, (1) and (2) if we set y = v » 0, but decay in absolute 

value if y 0, V # 0, The decay of |h I need not be mono- 

c 

tonic, but often is; the decay of E is monotonic, 

THE COMPUTATION 

The simulation results described here consist of fovir runs, 
designated A, B, C, and D, All four were obtained using a 
(128) code; that is, wave vectors with integer components 
ranging from k^^^ = 1 to k ^^^ = 6k and their negative counter- 
parts are included. Initial nonzero |B(k,0)|^ and |y(k,0) 

~ ~ 2 ~ 

were confined to a ring in k space between k -9 and k =25* 
For runs A, B and C there was twice as much magnetic energy 
as kinetic energy in each mode, but for run D, kinetic and 
magnetic energies were equal. In all runs the spectrum was 
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independent of k within the ring and phases of magnetic modes 
were randomly chosen. The amount of cross helicity initially 
present is determined by selecting pheise euigles for kinetic 
Fourier amplitudes to make a specified angle 6 with respect 
to the magnetic amplitude at the same This controls both 
the modal cross helicity, H^(k) “ Re y*(k)‘B(^) and the total 
cross helicity. Runs A and D are initially highly aligned, 
with cos”^(.95) * 9 while runs B and C are only moderately 
correlated with cos~^(.3) = 9. In runs A, B and D, the param- 
eters y * V = l/400, and in run C, y = v = 1/1000. Times are 
measured in units of Alfven transit tim^s of unit distance, 
with timesteps of 1/256 such units. Other runs, not described 
here, have shown results similar to these four. 


PHYSICAL RESULTS 


IVpical behavior is exhibited in Fig. 1: a monotonic increase 

of |h /e|. If we start with H /E negative, it decreases with 
c c 


time . 

A simple explanation is suggested by writing Eqs. 
(2) in the Elsasser variables Z“ = y ^ B. We assume 
moment that inteirest lies in the low-k region in the 
space, where the dissipative terms may be neglected, 
and (2) may be combined to give 

K 5* - 

z" » -z"^-V2“ - 7p 


(l) and 
for the 
Fourier 
Eqs. (1) 

(5) 

(6) 


The pressTire p is, as usual, obtained by taking the divergence 

4> 

of either Eq. (5) or Eq. (6) and using V*Z“ = 0 to yield a 
Poisson equation. 

Now 8upix>se that in some mean square sense, Z* is large 
con^ared to Z~ ■ : IZ""! » |2~|. It follows that Z* is nearly 
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time independent, since its time derivative is of 0 (Z~): 

(5a) 

[7p is of the same order as g *72^]. Hovever, the evolution 
of Z is not slow, £ind is given to lowest significant order 

Z" 5 -2‘"(t=0)-7z‘ - 7p, (6a) 

where 2^(0) ia not negligible. Thus the fractional rate 
of change |2*| ^|9Z*/9t| is small, but |z”|”^| 9Z~/9t | is not. 

No spectral transfer in the Z* field is implied by Eq. (5a). 
Its spectrvun is frozen at this order; if it is initially con- 
fined to wave numbers where dissipation is negligible, it will 
remain there. However, at every time step, Eq. (6a) will 
transfer the 2 spectrum to all additive combinations of wave 
numbers that can be made up of those present in Z*(0) and Z” 
at the previous time step. Clearly, the transfer of the Z” 
spectrum to higher wave numbers will continue at a fractional 
rate that does not go to zero eis Z” becomes small. Z” will be 
dissipated as it reaches the hi^er k values, thus enhancing 
the inequality jz*| » |z”|, (A somewhat different result 
was obtained by Dobrowolyny et al, 198O). 

We believe this mechanism, by which the majority species 
cannibalizes the minority one by sending it to high wave num- 
bers to be dissipated, is the essential mechanism involved in 
dynamic alignment. [A pure Z* field is aligned and a pure Z" 
field is anti-aligned, j A more complicated dynamical calcu- 
lation is required to make the effect precise, particularly 
for 0(1) initial ratios of l2 ^|/|z”l* 

Figure 2 shows the normalized cross helicity spectnan 
2H^(lj)/E(ij) vs. k at an early stage of the evolution of run A, 
t as 0.19< The initial vcQ.ue of this quantity was 0.895, necu* 
to its maximum possible value of 1.0, in the range of nonzero 
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excitations, 3 £ k £ 5. Figure 2 shows that the transfer of 
excitation to higher k's than those initially present consists 
almost entirely of negative cross helicity contributions, 
corresponding to hi^Iy anticorrelated small-scale magnetic 
and velocity field fluctuations. This typical behavior shows 
the transfer of Z” to higher k \rtiere nothing was initially, 
and being dissipated there. 

Ihe picture seems, however, not to be as simple as the pre- 
ceding argument would indicate. A competing process, the 

"selective decay" of energy relative to mean square vector 

2 2 

potential A = Sj^|g(fe,t)| /k edso occurs under some circum- 
stances, and at times competes with the increase of |h^/E| 
illustrated in Figs. 1 and 2. A is also an ideed invariant, 
and in several other computations (Matthaexis and Montgomeiy, 

1980), the ratio E/A has been observed to decay monotonicsdly, 

2 

to close to its minimvm value of k . =1. 

min 

Ihe alignment process and selective decay cannot simulta- 
neously proceed to their limiting states, since the limits 

I2H /E| -►I and E/A k , are incompatible. Figure 3 shows 
' c min 

that for initially hi^ly aligned runs, A and D, the ratio 
(1 - 2H^/E)/(1 - A/E) decreases indicating that the alignment 
process dominates (note that = l). However, the ratio 

increases for the less aligned runs B and C. \t this writing, 
we do not know of any general theoretical considerations which 
determine which process "wins" in the long time high Reynolds 
nusiber limits, for arbltraiy initial conditions. 

We may conjecture a (complicated or simple) boundary in 
the phase space of initial conditions on opposite sides of 
which the tendency is toward alignment and toward a minimum 
energy state. Similar uncertainties exist with regard to the 
three dimensional case, where magnetic helicity takes over 
the role of A, and the final selectively-decayed state is the 
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forces free "Taylor state.* We have not succeeded, at this 
point, in identifying the boundary. 
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nOURE CAPTIONS 

Fig. !♦ H^/E for runs A and E vs. time in Alf-zln transit 

times of \mit distance. The increeise of K /E with time shows 

c 

that dynamic alignment of v and B is proceeding. 

Fig. 2. Normalized cross-helicity snectrum 2H (k)/E(k) for 

c 

run A at 5 = .19. H^(k) euad E(k) are modal spectra, i.e. 

average values over all k with length k. Initially H (k) = 

^ c 

.895 for 3 £ k £ 5 and is zero elsewhere. Spectrad transfer 
to hi^er k is shown to be almost entirely negative cross- 
helicity excitation (that is, the Z” field). 

3 » (1 - 2H /E)/(1 - A/E) for runs A, B, C, D. Increase 

c 

of this ratio indicates dominance of selective decay; de- 
crease signals dynamic alignment. The suggestion is that 
initially highly aligned initial conditions leads to ultimate 
dominance of alignment, but no threshold has been established. 


